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n  index  of  refraction 
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q  electronic  charge 
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v  fluid  velocity  relative  to  the  standing  shock  (lab) 

W  atomic  weight  of  argon  atom 

W  average  atomic  weight  of  the  mixture 

x  degree  of  ionization 

Greek  Letters 

OC  polarizability 

7  ratio  of  specific  heats 

A(  )  change  in  (  ) 

5  shock  stand-off  distance 

9  argon  ionization  temperature 

X  observing  wave  length;  also  X  =  4500A 

X'  4030A 

M-  reduced  mass  of  ion- electron  pair 

vii 


NOTATION  (Continued) 

p  mass  density 
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I .   INTRODUCTION 

Treatment  of  high  temperature  shock  layer  flows  has  been  of  interest 

from  the  point  of  view  of  hypersonic  and  re-entry  heat  transfer.   Fay 

r  i*) 
and  Riddell  LI J   developed  a  correlation  formula  for  convective  stag- 
nation point  heat  transfer  rate  in  a  dissociating  diatomic  gas,  in 
conjunction  with  an  experimental  study  by  Rose  and  Stark  [2].   Several 
investigators  have  applied  a  modification  of  this  formula  to  stagnation 
point  heat  transfer  in  partially  ionized  argon.  Reilly  [3]>  operating 
in  a  regime  where  radiative  heat  transfer  was  considered  negligible, 
measured  convective  heat  transfer  rates  slightly  higher  than  predicted 
by  the  Fay-Riddell  equation  applied  to  an  ideal  gas.  Rutowski  [h]   also 
determined  that  convective  heat  transfer  was  only  moderately  increased 
by  real  gas  effects  in  the  partially  ionized  argon.  Rutowski  and  Ber- 
shader  [5]  found  that,  for  a  shock  Mach  number  of  twelve,  radiative  and 
convective  heat  transfer  rates  were  about  equal,  while  the  convective 
rate  alone  agreed  with  the  Fay-Riddell  ideal  gas  calculation.  They  were 
also  able  to  correlate  the  measured  radiative  heat  transfer  to  the  cal- 
culated electron  density,  using  the  Kramers  classical  theory  to  relate 
radiation  to  the  square  of  electron  density. 

The  importance  of  electrons  in  both  radiative  and  convective  ther- 
mal transport  leads  to  interest  in  electron  density  measurement.  The 
present  study  deals  with  a  measurement  of  electron  density  near  the 
stagnation  line  of  a  two-dimensional  shock-heated  argon  flow  around  a 
cylinder-plate  (Fig.  l).  The  purpose  of  the  work  was  to  make  a  quanti- 
tative comparison  between  measured  results  obtained  by  optical  inter- 
ferometry  and  the  predictions  of  equilibrium  flow  theory,  using  the 
Rankine-Hugoniot  and  Saha  equations.  A  noble  gas  is  desirable  as  the 
working  fluid,  since  ionization  can  be  studied  without  the  complications 
of  dissociation.   In  particular,  argon  was  chosen  because  of  its  moderate 
ionization  potential  and  relatively  high  molecular  weight,  which  enabled 
significant  ionization  to  be  obtained  behind  the  standing  shock. 




'Numbers  in  brackets  subsequently  designate  References  collected 
at  the  end  of  this  paper. 


The  technique  of  optical  interferometry  depends  upon  changes  in 
index  of  refraction  of  the  test  gas,  and  has  long  been  used  to  map 
densities  in  non- ionized  flows.   But  Alpher  and  White  [6]  demonstrated 
that  interferometry  at  optical  wave  lengths  could  be  used  to  measure 
electron  densities  in  a  plasma.   Since  the  refractivity  of  atoms  is 
relatively  insensitive  to  wave  length,  while  the  refractivity  of  elec- 
trons is  highly  dispersive,  measurement  of  fringe  shifts  at  two  wave 
lengths  enables  the  electron  contribution  to  be  separated.   This  yields 
a  direct  measurement  of  electron  density,  independent  of  the  state  of 
equilibrium  of  the  plasma.   Klein  [7]  applied  this  method  to  two-dimen- 
sional flow  around  a  model  in  an  electromagnetically-driven  shock  tube, 
but  the  flow  entering  the  standing  shock  was  highly  non-uniform.   Wong 
[8]  has  utilized  this  two  wave  length  technique  to  study  the  relaxation 
zone  behind  the  moving  shock  in  argon.   His  work,  using  the  same  combus- 
tion driven  tube  as  the  present  study,  reinforces  the  conclusion  drawn 
from  Figs.  2  and  3  that  flow  entering  the  standing  shock  is  quite  uni- 
form. 

The  present  study  is  part  of  a  long  range  program  at  Stanford 
dealing  with  the  dynamics  and  properties  of  shock-heated  plasmas.   In 
particular,  the  present  problem  has  direct  application  to  ionization 
phenomena  in  the  close  environment  of  hypersonic  vehicles,  such  as 
communications  blackout. 


II.      EQUILIBRIUM  FLOW  THEORY 

A.  Introduction 

In  order  to  estimate  the  stagnation  region  electron  density,  it 
is  necessary  to  calculate  the  state  properties  of  argon  at  the  stagnation 
line.   Beginning  with  the  state  of  the  gas  in  the  expansion  chamber  prior 
to  the  arrival  of  the  incident  shock  wave  (Fig.  1,  region  l),  the  one- 
dimensional  state  calculations  are  made  across  the  incident  and  standing 
shocks,  assuming  adiabatic  flow  and  the  attainment  of  local  thermodynamic 
equilibrium  after  each  shock.   The  equilibrium  assumption  is  justified  for 
region  2  by  the  work  of  Wong  [8],  as  well  as  the  interferograms  in  Figs. 
2  and  3  of  the  present  work.   The  computation  is  completed  by  assuming 
that  the  gas  flows  isentropically  and  incompressibly  from  the  standing 
shock  to  the  stagnation  line. 

B.  Gas  Properties  Across  Shock  Waves 

Properties  across  both  shock  waves  are  calculated  using  the  Rankine- 
Hugoniot  equations  for  a  partially  ionized,  monatomic  gas  in  conjunction 
with  the  Saha  equation.   The  validity  of  the  Saha  equation  for  the 
equilibrium  argon  plasma  has  been  demonstrated  by  Wong  [8]  and  Cole  [12]. 
Results  of  a  computer  program  obtained  by  Wong  were  used  to  obtain  the 
state  of  the  gas  in  region  2.  The  equations  used  across  the  standing 
shock  follow. 

The  degree  of  ionization,   x  ,  is  defined  by 

N 
e 
x  = 


N  +  N 

A         A 


where  N,  ,   N    and  N   are  number  densities  of  atoms,  ions,  and 
A     A+ 

electrons  in  any  region  of  flow.   The  Saha  equation,  applied  to  the 
atom- ion-electron  gas,  is 

BA(T)ph3  e*/kT 

+  1 


2V 


(T)(27Tu)3/2(kT)5/2 


where  p  and  T  are  the  pressure  and  temperature,   h   is  Planck' s 
constant,   x   is  the  ionization  potential,  and  u  is  the  reduced  mass 
of  the  ion-electron  pair.   In  order  to  make  use  of  the  Saha  equation, 
it  is  necessary  to  approximate 


A 


and  B   ,  the  partition  functions 
A+ 


of  the  atom  and  ion.   Since  the  neutral  atom  of  argon  has  about  60 
excited  states  near  the  ionization  potential,  the  partition  function  is 
approximately 

162,500 
BA(T)  =  1  +  60e     T       , 

where  162,500  K.   is  the  ionization  temperature.   Since  the  temperatures 

involved  are  less  than  15,000°K. ,  it  is  evident  that  the  second  term  can 

be  neglected.   The  argon  ion  has  a  ground  state  of  multiplicity  four,  and 

a  low  level  excited  state  of  multiplicity  two.   Neglecting  higher  levels, 

the  partition  function  of  the  ion  is 

2060 

B   (T)  =  k  +   2e    T 
A 

Bond  [9]  estimates  that  neglect  of  these  higher  levels  and  the  second 
ionization  of  the  atom  will  not  introduce  appreciable  error  below 
15 > 000  K.   Using  these  approximations,  the  Saha  equation  becomes 

1-1/2 


x  = 


ph; 


182,100 
T 


/     2060  \ 
(27T  u)3/2(kT)5/2  k\2+e       T   j 


+  1 


The  equation  of  state  of  an  ionized  monatomic  gas  in  thermodynamic 
equilibrium  is 

p  =  (l+x)  pRT  =  pRT    , 


where  R  is  the  average  gas  constant  per  unit  mass  of  the  mixture,  and 
the  masses  of  the  ion  and  atom  are  considered  equal. 

The  one-dimensional  equations  of  conservation  of  mass,  momentum, 
and  energy  for  flow  across  a  normal  shock  are  unchanged  by  ionization: 


P2v2  =  p3v3 

P2V2  +  P2  =  P3V3  +  P3 

\  4  +  h2  =  \  v3  +  h3    ' 

The  velocities   (v)  are  given  in  the  laboratory  system,  i.e.  a  system 
fixed  with  respect  to  the  standing  shock.   These  equations  may  be 
combined  to  obtain 

vh2^|-^)(v2)   •  (n-1J 

If  excitation  energy  is  neglected,  the  specific  enthalpy  is  given  by 

h  =  eA  +  e .  +  ^    , 
tip7 

where  e   is  the  internal  energy  per  unit  mass  due  to  translation,  and 
e.   is  the  specific  ionization  energy.   They  are  given  by 


1 


3  —  3  R*  3  R*    /        \  3      / 

6t       2ET  =  2-     T       2F   (1+X)T  =  2  R(1+X)T        » 
W 


and 


x(N     +  N   )x       N  k9       N  k6 

e.    =  £ ^—  =  -£—  =     e ,  T    ,    v    =  xR0        , 

l  p  p  mjf/x) 

where  R*  is  the  gas  constant  per  mole,   W  is  the  atomic  weight  of 
argon,   W  is  the  average  mixture  molecular  weight,  9     is  the  ioniza- 
tion temperature,   and  m   is  the  mass  of  the  argon  atom.   Finally,  the 
specific  enthalpy  is 

h  =  |  (l+x)RT  +  xR0 
Substituting  this  in  Eq.  (il-l),  the  Rankine-Hugoniot  equation  becomes 


k-±  +  1 


p3  +  ,.2!£!s)£ 

l+x„  I  T. 


(H-2) 


A  third  equation,    obtained  by  applying  the  equation  of  state  on  each 
side  of  the  shock  wave,  is 


1+x3/p2  p3 


(H-3) 


A  fourth  independent  relation  is  found  by  eliminating  v   from  the 
continuity  and  momentum  equations  to  obtain 

2   P2   P3    (P3/P2^  "  -1 
V2  "  P2  V2        (P^/P2)  -  1 

The  densities  are  eliminated  from  the  preceding  equation  by  substitution 
from  Eq.  (II- 2)  and  the  equation  of  state.   The  result  is 


\h  +  1\  h-S 


v2  =    (l+x2)RT2 


3 


1+2 


X3~X2\  9 
1   +xj  T2 


(II-M 


To  obtain  this  equation  in  terms  of  Mach  number,  the  speed  of  sound  in 
region  2,   a, 

c 

expressed  as 


region  2,   a  ,   is  required.   For  a  monatomic,  ionized  gas,   a   is 


a  = 
2 


^~2-   J!  R(1+X2)T; 


Dividing  Eq.  (II-4)  by  a   ,  and  solving  for  the  pressure  ratio  Po/Po  > 
the  result  is 


'3       1 


zf  =  tj  <  3   5NU  - 
2 


P   Bv  ■  ""2 


h  M2,2  ,  32  (X3-X2)  0M2 

(1-M2}  +T5-n^g—2 


,~1 


1/2 


where  M   is  the  Mach  number  entering  the  standing  shock.   This 


equation,  along  with  the  Saha  equation,  and  Eqs.  (II-2)  and  (II-3), 
form  a  determinate  system  for  the  four  unknowns  p  ,   p   ,   T   and 
x   .   The  one-dimensional  conservation  equations  used  in  this  derivation 
are  valid  when  applied  to  equilibrium  regions  on  each  side  of  a  thin 
normal  shock.   Since  the  thickness  of  the  non- equilibrium  region  in 
Fig.  3  is  about  half  of  the  shock  stand-off  distance,  the  assumption  of 
one-dimensional  flow  through  the  shock  may  be  questionable. 

C.   Stagnation  Line  Gas  Properties 

If  the  flow  from  region  3  to  the  stagnation  line  is  assumed  adia- 
batic  and  non-heat  conducting,  then  the  total  enthalpy  is  constant  along 
any  equilibrium  streamline.   In  particular,  along  the  central  stream- 
line, the  enthalpy  is  given  by 

hs  ■  h3  +  \A         ■ 

The  integrated  one-dimensional  conservation  equations  used  across  the 
shock  are  inappropriate  as  the  flow  becomes  two-dimensional  approaching 
the  stagnation  line.   But  Euler' s  equation  along  the  central  stream 
tube  can  be  written  as 

vdv  +  £P  =  0 
P 

where  neglect  of  viscous  stresses  is  justified  by  the  relatively  high 
densities,  hence  Reynolds  numbers,  involved.   The  assumptions  of  adia- 
batic,  non-conducting,  equilibrium,  non- viscous  flow  imply  the  assump- 
tion of  isentropic  flow.   Integrating  the  previous  equation,  we  obtain 

pvi: 

pvdv  -  \   dp  =  ^  +   pg  -  p3  -  0 

If  the  flow  is  also  assumed  to  be  incompressible,  then  p=p  =const.  , 
and  the  result  is  simply  Bernoulli's  equation: 

Ps  =  p3  +  \   P3V3    • 


The  assumption  of  incompressible  flow  is  hardly  tenable  for  the  hot, 
ionized  gas,  even  at  the  low  Mach  number  involved.   But,  since  Rutowski 
[k]    estimates  the  density  change  to  be  small,  Bernoulli's  equation  will 
be  used  in  these  calculations. 

With  stagnation  enthalpy  and  pressure  known,  the  previously  derived 

equations  of  specific  enthalpy  and  state,  and  the  Saha  equation  are 

applied  at  the  stagnation  point  to  solve  for  p  ,   x  ,   and  T 

s     s  s 

These  equations  are 


h  =  ^  (l+x  )RT  +xR0 
s   2     s   s    s 


P_  = 


s  '  (l+x  )RT 
s   s 


and 


x  = 
s 


182, 100 

v.3    Ts 
p  h  e   b 

(2  tt  m)3^  )5/S  M2+e-206°/T=) 


+  1 


-1/2 


Results  of  the  iterative  solution  were  used  to  calculate  N    from  the 

e 
s 

relation 


P  * 

s  s 

s     A 


R 


III.   APPLICATION  OF  INTERFEROMETRY  TO  ELECTRON  DENSITY  MEASUREMENT 
The  Mach-Zehnder  Interferometer 


A  schematic  of  the  interferometer  is  shown  in  Fig.  9°  Monochromatic 
light  from  the  source  is  collimated  before  entering  the  half- silvered 
beam  splitter  S   .   Here  the  light  is  split  into  two  beams,  which  travel 
via  the  fully  silvered  mirrors  M   or  M   to  be  recombined  at  the 
second  beam  splitter  S   .   The  light  is  focused  by  the  lens  L  ,  and 
the  camera  is  placed  at  the  focal  plane  of  the  test  section. 

If  the  mirrors  were  all  perfectly  parallel,  and  the  optical  path 
lengths  of  the  two  beams  exactly  equal,  one  would  expect  reinforcement 
of  the  in- phase  beams,  resulting  in  a  uniformly  bright  image  at  the 
camera.   Now,  if  one  of  the  mirrors  is  rotated  slightly,  a  uniformly 
varying  difference  in  optical  path  is  introduced,  resulting  in  alter- 
nating reinforcement  and  destruction,  i.e.,  a  fringe  pattern.   If  now  a 
plasma  is  formed  in  the  test  section,  the  index  of  refraction,  and  hence 
the  optical  path  length,  of  one  beam  is  changed,  resulting  in  a  fringe 
shift  given  by 

A  nL   A  (n-l)L 
S  =  ~T~  ~     X 

where  A  (n-l)   is  the  change  in  refraetivity,   L  is  the  test  section 
width,  and  X  is  the  monochromatic  observing  wave  length.   The  com- 
pensating chamber  is  introduced  to  compensate  for  the  windows  of  the 
test  section  of  the  shock  tube. 

Since  the  primary  quantity  desired  is  electron  density  at  the 
stagnation  line,  the  problem  is  to  relate  the  observed  fringe  shift 
(from  ahead  of  the  standing  shock,  2,  to  the  stagnation  line,   s  )  to 
the  change  in  electron  density.   Conditions  behind  the  moving  shock,   2  , 
are  taken  from  calculations  by  Wong  [8],  which  were  verified  experi- 
mentally by  him. 


B.    Development  of  Fringe  Shift  Equations 

In  order  to  relate  the  observed  fringe  shifts  to  electron  density, 
the  refractivities  of  the  neutrals,  ions,  and  electrons  must  "be  related 
to  thermodynamic  quantities. 

The  refractivity  of  the  neutrals  is  obtained  from  the  polarizabi- 
lity  of  the  atom,  a     ,  by  the  equation 

2tt  a 
nA  -  1  .  ar  Q :* A  =  — -   (H  m  )  =  ^  I     , 

A 

where  It  -  (2tt  a  /m   )   is  the  Gladstone-Dale  constant,  which  varies 

only  slightly  with  wavelength.  The  atom  number  density  is  related  to  the 
degree  of  ionization  by 


N.  =  £-  (1-x) 
A   mA 

Therefore,  the  change  in  refractivity  of  the  neutrals  between  any  two 
regions  is 

A  (nA-l)  =  K^A  [p(l-x)] 

Similarly,  the  ion  refractivity  is  given  by 

27T  a  +  a  + 

n    ,    -   1  =  27T  a     N       =  — •   —   •   a.m.N    ,    =  It    -—     m.N 

a+  »+   *+         OL.  m.  A  A    .+  A  a.  A    .+ 

A  AAAA  A  AA 

The  number  density  of  the  ions  is 

N   =  N  =  &        . 
A+    6   mA 

The  change  in  refractivity  of  the  ions  is 

V 

A  <nA+  -  1)  ■  \  ej"  A  (px)    • 


10 


The  refractivity  of  the  electrons  is  given  by  Wong  [8]  to  be 
2 


u.;T 


(ne-D  =  -  g   ~ 
where  u>  ,   the  plasma  frequency,  is  given  by 


n        Kir   q  N    ,    2 
e   47T  q 


WP=   m 


m  mA 
e  A 


(px) 


10  is  the  observing  light  frequency,  and  where  the  assumption  has  been 
made  that  the  electron-atom  collision  frequency  is  much  less  than  the 
observing  frequency.   This  assumption  is  valid  for  conditions  of  the 
present  study.   The  change  in  refractivity  due  to  electrons  is 


2.2 


A  (n  -1)  =   "  qX    A  (px)  =  -  6.77  x  10 VA  (px) 

2m  m.TTC 
e  A 

where  X  is  in  centimeters,  and  A  (px)   signifies  the  difference  in 
the  product  of  density  and  degree  of  ionization  between  any  two  regions, 

Assuming  that  the  refractivity  changes  for  the  various  particles 
are  independent,  the  total  change  in  refractivity  is 

A  (n-1)  =  A  (n  -1)  +  A  (n   -  l)  +  A  (n  -l)   . 

A         A  e 

Recalling  that  the  total  fringe  shift  across,  for  example,  the  standing 
shock  is  given  by 


S  =  J  A  (n-1) 

this  fringe  shift  is  expressed  in  terms  of  thermodynamic  quantities  by 
r 


S^<r 


V 

P3(l-x3)-p2(l-x2)+  —  (P3*3-P2x2) 


k- 


•  6.77  x  108X(p3x3-p2x2); 


(III-D 
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The  ratio  of  polarizabilities  of  the  argon  ion  and  atom  is  taken  as 

a 
A+ 

cT"  =  °'67 

A 

and  the  test  section  width  for  the  present  experiment  was  L=5.0  cm. 
Inserting  these  values  and  combining  terms,  the  total  fringe  shift  for 
X  =  4500  A  is 

S  =  1.78  x  io4[p3-p2-8.86(p3x3  -  p2x2)]    ,  (II1-2) 

where  K.  ,_   =  .1592  cm  /gm  has  been  used.   Similarly,  for  X'=4030  A  , 
K.  _  =  .1609  cm  /gm  ,  and  the  fringe  shift  is 

S'  =  1.993  x  104[p3-p2-7.13(p3x3  -  p2x2)]   .  (III-3) 

Combining  the  last  two  equations,  the  electron  density  in  region  3  of 
the  argon  plasma  is 

N   =  N   +  4.33  X  101T(S'  -  1.1275)     ,  (UI-^) 

3    e2 

where  X  =  ^500  A  ,   X1  =  4030  A  ,   and  S  and  S'   are  the  measured 

fringe  shifts  at  these  wave  lengths.  With  N    now  determined,  either 

63 

of  the  fringe  shift  equations,  (III-2)  or  (III-3),  can  be  used  to  solve 

for  p   .   Choosing  the  equation  for  X  =   ^500  A  ,  the  result  is 

p0  =  p  +  5.63  x  10"5  S  +  5-86  x  10"22(N   -  N   I 
3    2  [   e3    ej 

The  degree  of  ionization  can  now  be  determined  from 

N  m. 
e  A 

x_  =  — =2 — 

3     p3 

The  last  three  equations  are  sufficient  to  determine  electron  density 
and  degree  of  ionization  at  any  point  in  the  plasma  from  observed  fringe 


12 


shifts  between  a  known  region  and  the  region  of  interest.  The  fringe 
shift  equations  (III-2)  and  (III-3)  may  be  used  to  predict  fringe  shifts 
from  calculated  plasma  conditions,  as  in  Fig.  h. 


13 


IV.   DESCRIPTION  OF  EXPERIMENT 

The  experiment  was  performed  using  the  Stanford  aerophysics 
laboratory  shock  tube.   Figure  8  is  a  schematic  of  the  experimental 
arrangement.  The  25  foot  long  driven  section  is  constructed  of  extruded 
aluminum,  with  a  two  inch  square  inside  cross  section.   The  combustion 
driver  uses  a  stoichiometric  mixture  of  hydrogen  and  oxygen  diluted  with 
helium,  which  is  ignited  by  a  battery  of  spark  plugs.  The  driver  and 

driven  sections  are  separated  by  a  scribed  aluminum  diaphragm.   Pumping 

-5 
of  the  driven  section,  to  about  10   mm.  Hg  before  filling,  is  accom- 
plished through  a  station  at  each  end  of  the  section.  The  test  section, 
located  20  feet  from  the  diaphragm,  supports  the  two-dimensional  cylinder- 
plate  model.  A  pair  of  optical  windows  permits  viewing  of  flow  about  the 
model,  and  three  barium  titanate  piezo-electric  transducers  just  upstream 
of  the  model  are  used  to  measure  the  incident  shock  speed.   The  dump 
tank  at  the  downstream  end  is  separated  from  the  driven  section  by  a  thin 
aluminum  diaphragm  which  breaks  on  arrival  of  the  incident  shock. 

The  Mach-Zehnder  interferometer  (Fig.  9)>  located  at  the  test 
section,  is  on  loan  through  the  courtesy  of  Lockheed  Missiles  and  Space 
Company.   The  light  source  was  obtained  by  discharging  l6  kilovolts 
from  a  7-5  microfarad  capacitor  across  magnesium  electrodes  spaced  about 
l/h   inch  apart.  Although  no  direct  measurement  of  the  spark  duration 
was  made,  oscilloscope  pick-up  of  the  capacitor  "ringing"  indicated 
that  the  spark  was  probably  near  maximum  intensity  during  the  entire 
exposure  time  of  five  microseconds.  A  blast  shutter  and  pinhole  were 
originally  placed  at  the  focal  point  of  the  light  source  in  an  attempt 
to  reduce  film  exposure  by  plasma  radiation.   However,  Fig.  10  shows 
that,  even  at  relatively  low  Mach  number,  radiation  reaching  the  film 
was  intense  enough  to  wash  out  completely  the  fringes  through  the 
5890A  filter. 

Since  it  had  become  evident  that  the  mechanical  blast  shutter  was 
not  fast  enough  for  this  application,  a  Kerr  cell  was  obtained,  again 
through  the  courtesy  of  Lockheed.  This  electronic  shutter,  manufactured 
by  the  Kappa  Scientific  Corporation  of  Sierra  Madre,  California,  was 
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equipped  with  a  pulse- forming  network  that  permitted  a  very  reproducable 
shutter-open  time  of  five  microseconds.   For  the  present  steady  flow 
application,  this  exposure  time  proved  satisfactory.   Unfortunately  for 
this  application,  it  was  designed  for  use  in  the  ultraviolet,  with 
acetronitrile  as  the  working  fluid.  The  Kerr  cell  was  triggered  to  open 
for  five  microseconds  during  peak  spark  intensity,  and  the  shutter-open 
period  was  always  between  12  and  17  microseconds  after  passage  of  the 
incident  shock  by  the  model. 

The  flow  field  was  split  by  interference  filters  passing  4500  and 
5890  angstroms  located  just  ahead  of  the  Polaroid  camera.   Simultaneous 
measurements  of  near  stagnation-line  fringe  shifts  through  the  two  fil- 
ters were  sufficient  to  determine  the  electron  density.  However, 
variation  in  intensity  of  the  dark  fringes  passed  by  the  589OA  filter 
(see  Figs.  2,3)  indicated  the  presence  of  a  harmonic.   Further,  compa- 
rison of  observed  fringe  shifts  with  Eqs.  (III-2)  and  (III-3)  showed 
that  the  "5890"  fringes  actually  corresponded  to  a  wave  length  lower 
than  U500A  .   Examination  of  a  spectrum  obtained  by  Wong  through  this 
second  order  589OA  filter,  also  using  a  magnesium  spark  source,  showed 
a  third  order  pass  band  at  4030A  .   It  was  concluded  that  peaking  of 
the  ultraviolet  Kerr  cell  at  about  4-000A  ,  as  well  as  the  low  intensity 
of  the  magnesium  spectrum  at  589OA  ,  resulted  in  a  fringe  pattern  domi- 
nated by  X'  =  4030A  .   The  existence  of  harmonics  in  the  fringe  pictures 
indicates  that  the  presence  of  some  589OA  might  cause  a  small  increase 
in  effective  wave  length  above  the  value  of  X1  =  4-030A  used  in  the 
calculations. 

An  attempt  was  made  to  obtain  total  stagnation  line  heat  transfer 
data  using  a  platinum  thick- film  calorimeter,  as  described  by  Rutowski 
[h] .      Initial  efforts  indicated  that  this  phase  of  the  program  would 
require  more  effort  than  originally  anticipated.   Because  of  practical 
limitations  of  time,  it  was  decided  to  concentrate  remaining  efforts  on 
the  interferometric  studies. 
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V.   DISCUSSION  OF  RESULTS 

Usable  fringe  shift  data  were  obtained  for  both  X  =  4500  A   and 
X'  =  ^030  A  on  four  runs.  These  were  made  at  incident  shock  Mach 
numbers  of  11.45,  11.8,  12.0,  and  12.8.   The  lowest  Mach  number  was 
obtained  with  an  initial  pressure   (p  )   of  10  mm.  Hg.  ,  while  the  others 
corresponded  to  P-.=5  mm.  All  tests  were  made  in  argon,  with  the  l/2 
inch  diameter  cylinder-plate  model  installed  in  the  test  section.   Graphs 
of  fringe  shift,  density  ratio,  electron  density,  and  degree  of  ioniza- 
tion (Figs.  4,5,6  and  7)   compare  predicted  and  measured  quantities  in 
region  3,  just  downstream  of  the  standing  shock.   Figure  7  also  shows 
measured  stagnation  line  degree  of  ionization  and  the  corresponding 
prediction  of  equilibrium  theory. 

Attainment  of  equilibrium  in  region  2,  the  flow  entering  the 
standing  shock,  was  found  by  Wong  [8]  to  occur  about  six  microseconds 
after  the  incident  shock  (lab  time)  for  the  lowest  Mach  number  of  the 
present  study,  while  the  prediction  of  Petschek  and  Byron  [10]  is 
somewhat  less.   Since  the  present  interferograms  were  taken  between  12 
and  17  microseconds  after  the  incident  shock,  it  is  reasonable  that  the 
flow  entering  the  standing  shock  was  in  equilibrium.   The  interferograms 
in  Figs.  2  and  3  show  relatively  straight  fringes  in  region  2,  supporting 
the  conclusion  of  equilibrium  there.  The  compensating  effects  of  inci- 
dent shock  wave  attenuation  and  radiative  cooling  on  the  flow  in  region 
2  were  discussed  by  Rutowski  and  Bershader  [5l»   Although  the  magnitude 
of  these  effects  has  not  been  determined,  their  effect  probably  was 
minimized  by  taking  all  the  interferograms  at  the  same  time  relative  to 
the  incident  shock.  Total  available  test  time  was  estimated  to  be  about 
100  microseconds  for  the  fastest  run. 

The  interferograms  at  the  two  lower  Mach  numbers  (Fig.  2a)  show  a 
sharp  upward  (positive)  fringe  shift  across  the  standing  shock  due  to 
change  in  refractivity  of  the  heavy  particles  [Eq„  (ill-l)].  These 
shifts  were  found  to  correspond  not  to  the  predictions  of  equilibrium 
theory  [Eq.  (III-2)],  but  to  the  estimate  obtained  by  assuming  degree  of 
ionization  frozen  at  the  value  in  region  2,  and  equilibration  of  the 
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translational  modes  only,   Predicted  quantities  based  on  these  assump- 
tions were  obtained  from  the  perfect  gas  shock  tables  of  Keenan  and 
Kaye  [ll],  for   7=5/3  >   and  the  results  are  compared  with  measured 
values  of  fringe  shift,  density  ratio,  and  degree  of  ionization  in 
Figs,  k,    5,   and  7°   Since  measurements  at  the  stagnation  line  showed 
only  a  small  increase  in  degree  of  ionization  from  that  in  region  3 
(Fig.  7)>  the  flow  field  from  the  standing  shock  to  the  stagnation  line 
is  concluded  to  be  definitely  out  of  chemical  equilibrium,  and  in  fact 
nearly  frozen.  Measured  values  of  electron  density  in  region  3>  which 
differed  only  moderately  from  those  measured  at  the  stagnation  line,  are 
shown  in  Fig.  6  to  be  much  less  than  the  corresponding  equilibrium 
predictions. 

The  interferograms  at  the  two  higher  Mach  numbers  (Fig.  2b)  show  a 
radical  change  in  the  flow  field.   The  shock  stand-off  distance  is  re- 
duced by  more  than  a  factor  of  two.  The  initial  upward  fringe  shift 
again  roughly  corresponds  to  the  frozen  prediction,  with  the  larger 
shift  at  the  lower  wave  length,  as  would  be  expected  from  Eq.  (ill-l), 
for  the  heavy  particles .   But  a  large  shift  in  the  opposite  direction, 
i.e.  due  to  electron  production  [Eq.  (III-l)]  is  also  observed.   Since 
this  shift  is  proportional  to  wave  length,  it  is  greater  for  X  =  4500  A  , 

The  flow  appears  to  approach  local  equilibrium  at  the  line  of 
maximum  downward  fringe  shift,  located  about  half  way  from  the  shock  to 
the  stagnation  line.  Since  the  fringe  shift  due  to  electrons  from 
region  2  to  3  is  somewhat  less  than  predicted  for  equilibrium  (Fig. h) t 
the  measured  values  of  electron  density  obtained  from  Eq.  (111-4)  only 
approach  the  equilibrium  estimate.   It  is  possible  that  application  of 
the  one-dimensional  conservation  equations  across  a  shock  extending  half 
way  to  the  stagnation  line  is  invalid,  which  could  explain  the  lower 
measured  equilibrium  densities.  A  better  theory  would  make  use  of  a 
one-dimensional  analysis  across  the  frozen  portion  of  the  shock,  and 
calculate  the  approach  to  equilibrium  by  a  two-dimensional  theory  utili- 
zing a  rate  equation.   However,  the  change  in  the  fringe  pattern  (Fig.  3) 
perpendicular  to  the  flow  direction  is  small  near  the  centerline  of 
flow,  indicating  that  the  flow  across  the  shock  is  nearly  one-dimensional 
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in  the  region  of  interest.  The  flow  from  region  3  to  the  stagnation 

line  causes  an  upward  shifting  of  fringes,  which  for  these  conditions 

can  only  be  caused  by  ion-electron  recombination  [Eq.  (III-2)].   This 

effect,  which  of  course  is  not  predicted  by  the  adiabatic  equilibrium 

analysis,  is  illustrated  by  a  small  decrease  in  degree  of  ionization 

from  region  3  to  the  stagnation  line  at  M  =12.0  ,  and  a  larger  decrease 

s 

at  the  higher  Mach  number  (Fig.  rj).     Although  the  Kerr  cell  effectively 
reduces  the  interferogram  exposure  from  plasma  radiation,  other  pictures, 
using  only  the  blast  shutter,  show  the  existence  of  strong  radiation  in 
this  Mach  number  range.   It  appears  to  be  concentrated  in  the  region 
behind  the  bow  shock,  and  is  probably  largely  due  to  radiative  recombi- 
nation.  It  is  concluded  that  the  flow  at  the  higher  Mach  numbers  does 
approach  equilibrium  momentarily,  but  immediately  begins  cooling  due  to 
radiative  losses.   Radiative  recombination  could  also  explain  the  lower 
than  expected  equilibrium  electron  densities  measured  in  region  3«   The 
sharp  change  in  stand-off  distance  over  a  very  small  Mach  number  range 
may  be  related  to  the  attainment  of  equilibrium  at  the  higher  speeds. 
Table  II  summarizes  these  observations  on  the  flow  interferograms. 

Although  no  rigorous  theory  has  yet  been  developed  for  convective 
stagnation  region  heat  transfer  in  a  partially  ionized  monatomic  gas, 
the  Fay-Riddell  correlation  equation  [l]  has  been  applied  to  this  situ- 
ation by  Reilly  [3]  and  Rutowski  and  Bershader  [5].   Predicted  heat 
transfer  rates  based  on  this  equation  depend  on  the  value  chosen  for 
the  Lewis  number,  but  do  not  vary  markedly  from  the  ideal  gas  prediction 
based  on  the  same  equation.   Both  Reilly  and  Rutowski  and  Bershader 
found  substantial  agreement  of  measured  convective  heat  transfer  rates 
with  estimates  based  on  the  equation  of  Fay  and  Riddell.   The  present 
study  suggests  that  the  flow  field  ahead  of  Reilly' s  model  was  probably 
nearly  frozen,  since  his  study  used  a  much  smaller  model  and  lower 
initial  pressures  than  the  present  work,  and  involved  Mach  numbers  of 
12  or  less.  A  similar  comparison  of  the  present  study  with  the  experi- 
ments of  Rutowski  and  Bershader,  involving  Mach  numbers  greater  than 
11.7*  initial  pressures  of  10  mm.  Hg.  ,  and  a  larger  model,  indicates 
that  equilibrium  was  probably  attained  in  the  stagnation  region  during 
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their  work.  These  considerations  suggest  that  the  state  of  equilibrium 
of  the  shoek-heated  argon  plasma  does  not  critically  affect  stagnation 
region  convective  heat  transfer.  This  result  is  not  entirely  unexpected, 
since  convective  heat  transfer  is  basically  a  boundary  layer  phenomenon. 
In  the  first  approximation,  it  should  be  relatively  insensitive  to  shock 
layer  rate  processes,  particularly  for  the  relatively  high  densities 
involved  in  the  present  study  (Table  I). 

Rutowski  and  Bershader  suggest  a  radiative  heat  transfer  model 
based  on  an  extension  of  Kramers  classical  theory  for  free-free 
(bremsstrahlung)  radiation  to  include  radiative  recombination  into  the 
upper  bound  states.  Assuming  thermodynamic  equilibrium  they  estimate 
radiative  heat  transfer  rate  to  be  proportional  to  the  square  of  elec- 
tron density.  Use  of  predicted  electron  densities  in  the  present  case 
would  have  resulted  in  grossly  overestimated  radiative  heat  transfer 
for  the  frozen  tests,  but  a  much  smaller  error  at  the  higher  Mach 
numbers,  where  equilibrium  was  attained. 

The  following  factors  were  considered  as  possible  contributors  to 
error  of  the  experimental  measurements: 

1.  Error  in  determination  of  Mach  number  by  the  pressure 
transducer-oscilloscope  system  is  estimated  to  be  about  1$. 

2.  Actual  initial  temperatures  varied  from  T  =297  K  to 

T  =30^  K  .   Since  all  theoretical  calculations  were  based  on 

T  =300  K  ,  a  small  scatter  in  the  (uncorrected)  measured 
values  would  be  expected. 

3«   The  work  of  Wong  [8]  indicates  that  the  presence  of  impurities 
could  cause  a  decrease  in  the  relaxation  time.   However,  this 
effect  was  predicted  primarily  for  very  low  initial  electron 
densities,  where  atom-atom  collisions  are  predominant.   Since 
all  regions  of  interest  in  the  present  study  involve  relatively 


was  about  10   mm,  Hg.   for  all  runs,  the  effect  of  impurities 


high  electron  densities,  and  the  tube  pressure  before  filling 
was  about  10~  mm,  Hg.   for  a. 
is  not  considered  significant. 
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k.        Fringe  shifts  could  be  read  to  an  accuracy  of  about  one-tenth 
of  a  fringe.  An  error  in  this  measurement  would  introduce 
only  a  small  discrepancy  between  predicted  and  measured  fringe 
shifts  at  a  single  known  wave  length  (Fig.  h,      X=4500A) . 
Therefore,  it  is  considered  that  qualitative  observations  from 
Fig.  k,    i.e.  separation  of  frozen  and  equilibrium  regimes,  are 
valid . 

5.   However,  the  measured  values  of  electron  density,  degree  of 
ionization,  and  density  ratio,  which  were  all  obtained  from 
measured  fringe  shifts  at  two  different  wave  lengths,  are 
subject  to  more  serious  error.  The  proximity  of  the  two  ob- 
serving wave  lengths   (X=4500A,   and  X'=ij-030A)   results  in 
sensitivity  of  Eq.  (III-4)  to  small  errors  in  reading  the 
fringe  shift,  especially  at  lower  Mach  numbers  where  the 
shift  due  to  electrons  is  small.   It  will  be  recalled  that 
original  plans  had  called  for  using  X'=5890A  as  the  second 
wave  length,  which  would  have  resulted  in  considerably  better 
accuracy.   There  is  an  additional  uncertainty  in  the  actual 
value  of  the  lower  wave  length  (X'=J+030A  +  ~100A  ,  -  OA)  . 
This  is  derived  from  the  possibility  that  presence  of  traces 
of  589OA  might  result  in  an  "effective"   X'  >  4030A  . 
Elimination  of  all  589OA  radiation  by  procuring  a  suitable 
additional  filter  was  precluded  by  time  considerations .   Use 
of  a  higher  "effective  wave  length  vice  X'=^030A  in  the 
calculations  would  have  resulted  in  a  fairly  uniform  increase 
in  the  measured  electron  density  for  all  runs.   The  estimated 
uncertainty  in  measured  electron  density  due  to  these  consi- 
derations is  included  in  Fig.  6.   It  should  be  noted  that  this 
interferometric  method  is  capable  of  considerably  greater 
accuracy  than  the  present  work  obtained.   However,  the  agree- 
ment of  density  ratios  (Fig.  5)  obtained  by  the  two-wave 
length  technique  with  the  appropriate  theory  (frozen  or  equi- 
librium) corresponds  to  the  conclusions  drawn  from  the  more 
accurate  single-wave  length  fringe  shift  predictions  (Fig.  k) . 
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VI.   CONCLUSIONS 

1.  Flow  entering  the  standing  shock  was  in  equilibrium  during  each 
test. 

2.  Stagnation  region  equilibrium  was  not  attained  during  the  two 
lower  Mach  number  tests.  The  degree  of  ionization  remained 
essentially  frozen  at  its  value  entering  the  standing  shock. 

3-   During  the  higher  speed  tests,  local  thermodynamic  equilibrium  was 
achieved  ahead  of  the  stagnation  point,,  The  thickness  of  the  re- 
laxation region  was  about  half  of  the  shock  stand-off  distance. 
The  degree  of  ionization  after  relaxation  was  less  than  the  equi- 
librium prediction,  and  actually  decreased  from  the  end  of  the 
relaxation  zone  to  the  stagnation  line.  The  lower  than  expected 
degree  of  ionization  in  both  regions  3  and  s  is  attributed  to 
radiative  recombination. 

k.        Since  radiative  heat  transfer  has  been  related  to  the  square  of 

electron  density,  the  present  work  indicates  somewhat  less  radia- 
tive heat  transfer  than  predicted  by  equilibrium  theory.   But  the 
observed  recombination  provides  some  justification  for  considering 
radiative  recombination  a  primary  component  of  stagnation  region 
thermal  transport. 
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X'=^030A 


X=1+500A 


Fig.    2a.      Flow  Interferogram     (M  =11.^5) 


X'=4030A 


X=4500A 


Fig.    2b.      Flow  Interferogram      (M  =11.8) 
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X'=ll-030A 


X=^500A 


Fig.    3a.      Flow  Interferogram     (M  =12.0) 

s 


X'=4030A 


\=4500A 


Fig.    3"b.      Flow  Interferogram      (M  =12.8) 
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FRINGE  SHIFTS  ACROSS  STANDING  SHOCK 
vs 
INCIDENT  SHOCK  MACH  NUMBER   (X  =  1+500A) 


Fig.  k.      Fringe  Shifts. 
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^— ^—  A  p  =  5  mm  Hg 

Q  p  =  10  mm  Hg 


DENSITY  RATIO  ACROSS  STANDING  SHOCK 

vs 

INCIDENT  SHOCK  MACH  NUMBER 

X  =  4500A 
X1  =  1+0  30A 


EQUILIBRIUM 
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M   (INCIDENT  SHOCK  MACH  NUMBER) 


Fig.  5.   Density  Ratio, 
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(cm"3) 1017 


A    p  =  5  mm  Hg 


□    p  =  10  ram  Hg 


ELECTRON  DENSITY  IN  REGION  3 

vs 

INCIDENT  SHOCK  MACH  NUMBER 

X  =  1+5  OOA 

X'    =  1+030A 
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Fig.    6.      Electron  Density. 
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Fig.  7«   Degree  of  Ionization. 
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X=1+500A 


Fig.    10.      Blast  Shutter  Interferogram. 
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